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METHODS FOR DETERMINING THE ORBITAL TRANS FE] 

REQUIRE}-!— — FOR NEAR-EARTH ORBITS 
By ‘Sanruel L. Miller and David J. Griffi" 

SUMMARY AND INTRODUCTION 

The Apollo C mission has several SPS test objectives that depend on 
Doth fixed and variable' burn durations for the rendezvous and deorbit 
maneuvers. In the process of defining the mission, the rendezvous 
maneuvers and the duration of SPS burns are subject to change, thus altering 
all planned maneuvers . This suggests a need for a method that quickly 
and accurately estimates the changed maneuver. This document provides 
both a graphical (plots) and an analytical method (Olivetti programs) for 
estimating the delta V requirements for the orbital transfer maneuvers . 

In addition -to the total delta V requirements , the components of. delta V 
in the external delta V system can be determined. 

With the exception of those for the rendezvous, the parametric data 
presented can be used to determine the orbital maneuver requirements. The 
use of the parametric data should thus reduce considerably the numbers 
of computer runs required to obtain .an exact solution. Both the Olivetti 
programs and the plotted data assume impulsive maneuvers; however, a good 
comparison of results with a finite burn can be obtained by placing 
the center of the burn arc at the impulsive point. 

• Use of both methods indicates that either gives comparatively 
accurate external delta V components. Thus, the mission planner has 
two independent sources of determining the orbital transfer requirements. 


Although this document was designed using the Apollo C mission 
profile, the data presented may be used for general mission planning. 
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ABBREVIATIONS 

Symbols 


a 

e 

h 

r 


.V . 



e 


semimajor axis 
eccentricity 
altitude 3 n.. mi. 
radius 

inertial velocity , fps 
total AV magnitude, fps 

components of AV in the external AV system, -fj 

pitch, i.e., the angle between the velocity vector and the* 
thrust vector measured in the orbital plane, positive up, deg 


y inertial flight-path angle, deg 

6 orbital central angle between perigee and vehicle position, deg 

y’ ' ' earth's gravitational constant, 1.U0T648 x 10 1 ' 6 ft 3 /sec z - 

iD yaw, i.e., the angle between the orbital plane ana me projec- 
tion of the thrust vector in a plane normal to the orbit plane, 
positive to the right', deg 


Subscripts 


a apogee 

E earth 

p perigee 

1 initial 


2 


desired 
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DISCUSSION OF METHODS 

i'r 11 ** 19 for specified 

a’ p 9 * ¥ > »* 

With the’ preceding parameters available, the targets for' a particular 
orosram^ & ° tained either *>7 the plotted data or hy the Olivetti .analytic 


The 


Figures 

using analytic equations. 


• Fi§ure 1 Presents the total velocity magnitude, AV^ , as a function 

2,ss .sir z sr ttoe - A ’ is ° inoiuaea in figure 1 *•■*** 

£ Z 

- P if Ur % 2 P + esent f J he inertial velocity and inertial flight-path ' 
ngle as a function of h^, h a> and 6 for various spacecraft altitudes 

ranging from 75 n. mi. to 265 n. mi., in increments of 10 n. mi. 

... Flgu ^ e 3 presents several nomograms giving the Z component of velo 
LLf® f function of present -inertial velocity and flight-path angle; 
pax^eterf 1V6S ^ X COmponent of velocity as a function of < the same 

Figure 5 presents the lateral .velocity component, AV y , as a functic 
of yaw and pitch for various AV^s and longitudinal velocity components ; 
AV y s. The magnitudes of AV^ presented are representative of those assc 
ciated with the SPS burns for the 205/101 Mission 


Olivetti Programs 

®ie following Olivetti programs may he used to obtain targets for 
a particular SPS burn without reference to plotted data. 
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Programs 1 and 2 calculate the velocity components in the orbital 
plane required for a particular transfer. , Program 1 uses both altitude 
and velocity vector inputs. Program. 2 uses only altitude inputs. 

Program. 3 calculates the lateral velocity component and the yaw ana 
pitch attitudes using the total velocity component and the in-plane veloc- 
ity components as input. 

Program k calculates the pitch orientation using the total velocity 
magnitude and the vertical velocity component as input. 

Progr am 5 uses the pitch attitude, the total velocity magnitude, 
and the in-plane velocity component -to calculate the lateral- velocity 
component and the yaw attitude. 

. These programs are presented in- appendices A througn Ji, respectively. 
The appendices provide instruction listings, input and output, equations 
used in each program, and the restrictions on the use of each program. 

Some. -equations given are not exact, but they dp yield an accurate approx- 
imation -.of the true value. 


Typical Example 

- - To -illustrate the use of the plotted, data and the Olivetti prpgrams , 

an example problem is presented. . Suppose the following is know: 

Present position Transfer ellipse and 


and 

ellipse 

. position desired 

h = 

175 mi. 

h = 175 n. mi. 

h 

= 114 n. mi. 

h . =121 n. mi 

p l 


P 2 

h 

= 228 n. mi. 

h - 277 n. mi 

a l 


a 2 

0 = 

97° 

o° < e 2 « i8o-° 


AV = 230 fps 
t 

V = 25 291 fps 
y = .898924° 

Find AV X , AV y , AV & , g, and ij>. 
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Graphical solution ,- The following method uses plotted, .data only: 


Step 1 - The inertial velocity and inertial flight-path angle on each 
ellipse is obtained from figure 2(k)T (These values are 
then used to determine AV^. and AVg,. ) 

V 1 = 25 296 fps V 2 = 25 395 fps 

' Y x = -91° Y 2 = 1-15° 

snep <z - Figures 3(c) and (d) are used to determine the velocity 

components in the Z direction, while figures Mb) and (c) 
give the velocity components in the X direction. 


V = 400 fps 
Z 1 


V = 510 fps 
• Z 2 



= »100 fps 


V T = 25 286 fps • 

X 1 

V Y = 25 391 fps 
A 2 

AY y = V v - = 105 fps 


Step 3 - The absolute value of pitch associated with the above AVg 
and AY, can be obtained from figure 1. 

~G 

| 6 | =29° 

Step 4 - From figure 5(a), the A Vy component and subsequently the 
yaw magnitude can be determined. 


• |ay y | = IT 3 fps 

I ip I = 59°- 

Analytic solution .- The following method uses Olivetti programs 1 
and 2 to calculate a solution. 

Programs 1 and 2 (appendices A and B) can be used to determine the 
AV„ and AVr, components directly from the data available - from the initial 

X z 

conditions . 
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Program 1 

Input: h. = IT 5 n. mi, * 

h • =- 121 n. mi. 
P 2 ‘ ' 

h ■ = 277 a- mi. 

a 

2 

A = 1 

Y = .89892it 0 
V 1 - 25 291 fps- 


Program 2 

- Input : h lJ5_.n ..".-ml:. 

h = llU n. mi 

P, . - . 


h- • = 228- n. mi 
a l 


h = 121 n. mi 
*2 

h = 277 n. mi 


Output: AV^. = 98.6026 fps Output: = 25 296 fps 

AV Z = -120.9121 fps | sin y^l = 0.0157381; 

V 2 = 25 392 fps 
* J sin Yg ] = 0.0203887 

AV X = 96 fps 
AF Z - -120 fps 


Using AV X and A7 Z 
can now determine AV Y> 


from program 1 as input, program 3 (appendix C) 
pitch and yaw magnitudes . 


Input; AV t •= 230 fps 

AV X = 98.6026 fps 
’ AV Z = -120.9121 fps 

Output: | j =^168.9906 fps 

|*| = 59.727^ 0 

|ej = 31.7158° - 

Comparison with numerical integration .- The example problem was 
solved using' numerical integration , and a comparison' with both plotted 
data and Olivetti programs is presented below. 
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'COMPARISION OF THE THREE METHODS 


[h = 175 

n. mi. ; h 

P 1 

= ll4 n. mi 

• 

# * 

228 n. ml. ; 

h = 121 n, 

. mi . ; h - 

277 el. mi 

p 2 

- 3 , _ 

2 




Methods 


Parameter 

Integration 

Figures 

Olivetti 

programs 

6, deg 

96.6 

95. 


\> fps 

25 291 

25 296 

25 291 

y ± > deg 

0.898924 

0.91 

0.898924 

e 2 , deg 

74.3 

73. 


Vg, fps . 

25 388 

25 395 

25 391 

Ygs deg 

1.17905 

1.15 

1.16818 

AV X , fps 

98 

105 

98.6 

AV y , fps 

166.346 

173 

168.99 

AV Z » fps 

-125 

-108 

-120.9 

AV. , fps 
X 

230 

230 

230 

i|» , deg 

59-4 

59. 

59.7 

8, deg 

32.8 

28. 

31-7 
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APPENDIX A 

I - TRANSFER VELOCITY COMPONENTS REQUIRED BASED ON 
BOTH ALTITUDE AND VELOCITY VECTOR INPUTS 


Macs**® 


*»* NOt 
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APPENDIX 'A 

. 

PROGRAM 1 - TRANSFER VELOCITY' COMPONENTS REQUIRED BASED ON 
BOTH ALTITUDE AND ■ VELOCITY VECTOR INPUTS 


A program is presented that calculates the velocity change necessary 
in the X and Z direction, to transfer impulsively from the present ellipse 

(h, V. , y. ) to the desired ellipse (h , h ) at some altitude h. 

11 a 2 Pg 

This program requires a decimal wheel setting of 7 with V being the entry 
point of the program. -The order of the input and output quantities is as 
indicated . 

Input: 'h, h , h_ , A (explained below), y," , V, 

. p 2 a 2 X 1 

The quantity A takes on the values of +1 or -1. If A = 1-, 
then the sign of y on the desired ellipse will be positive 
' (-i.. e . the vehicle will be on the. upleg. of the desired . 
ellipse after the impulsive maneuver). • If A =--l,'then the 
sign of y on the desired ellipse will be negative (i.e., the 
vehicle will be on the downleg of the desired ellipse after 
the impulsive maneuver). 

Output : AV X , AV Z 

Equations : 


r = h + r E 


r 

a 



+ r 


E 


r 

P 



+ r 


E 


a 


r 

a 


+ r 

E. 


2 
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AJ Z ~ Y 1 sin J .Y 1 - V 2 sin r 2 

Restrictions ; 

1. Neither the present ellipse nor ihe desired ellipse can he 
circular. 

2. This program should not he used when the flight-path angle of 
either ellipse'-exceeds 5° at the point of transfer, since a small angle 
approximation is used for both the sine -and 'cosine of y * 
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INSTRUCTION LISTING OF PROGRAM 1 


AV 

at 

•Rt 

rt 

R+ 

R+ 

V$ 

bt 

S 

i 

b+ 

B$ 

S 

4 - 

b+ 

bt 

S 

+ 

b+ 

at 

dt 

ct 

CT 


bt 

At 

At 

Ct 

bt 

Bt 

c- 

at 

'rt 

R- 

Rt 

Rt 

DX 

X 

A/ 

*4 

AX 

At 

bX 


bX 

Bt 

b$ 

A+ 

B- 

b$ 

bt 

At 

bt 

A/ 

CX 

s 

X 

bt 

a/ 

CX 

ct 

s 

t 

at 

Rt 

R- 


R* 

Rt 

Rt 

R- 

R* 

Rt 

TT 

D- 

AX 

Bt 

at 

at 

X 

•a/ 

Bt 

s 

X 

Bt 

X 

at 

RX 


rO 

R+ 

R- 

Rt 

Rt 

rX 

Rt 

RS 

DX 

ct 

Ct 

c| 

c- 

cX 

A<> 

bt 

cX 

st 

B- 

A<> 

V 
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appendix a 

PROGRAM 2 - TRANSFER VELOCITY COMPONENTS REQUIRED 
EASED ON ALTITUDE INPUTS 


IMPEDING 


PAGE BLAHK 130T. FUSSED. 
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APPENDIX B 

PROGRAM 2 - TRANSFER VELOCITY COMPONENTS REQUIRED 
BASED ON ALTITUDE INPUTS 


This program calculates the velocity change necessary in. -the X, and Z 
direction to transfer from the present ellipse (h , h *) to the- desired 

ellipse (h , h ) at an altitude h. A decimal -wheel setting of T-is 
P 2 a 2 

required for this program with V being the entry point of .the program. 

The order of the input and output quantities is as given. 

Input:, h, h , h , h , h 

P i ' a l P 2 a 2 

: Output : V lS .|sin yj,. V p » |sin yj» AV v ,. AV 

’ Equations : 


r = h + r. 


E 


r = h +' r_ 
p p -E 


r 

a 



+ r 


E 


a 


r 

a 


+ r 

R 

2 



V 



iin y| 




e 2 ) 

~rj 


CEDING 


bunk sot huso- 
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AV X W V 2 




V, 


sin 'y 



V, 


sin 


Restrictions : 

_ i 

1. Neither the present ellipse nor the desired ellipse can he 
circular'. " 


■ 2. The transfer from- the present to the desired ellipse must occur 
on the same side of: the line of apsides. 

3. AV -will have the correct sign, hut will have an. error in its 

X ... 

magnitude. The error in the magnitude of will increase as y increases 

in absolute value, since the cosine of y is not used in .the calculation. 

4. AV will have the correct magnitude and sign on the ascending 

Z 

side of the ellipse. On the descending side of the ellipse, the magnitude 

of AV„ will he correct hut its sign should he negated, 

z 

5. This program should not he used when the flight-path angle of 
either ellipse exceeds 5° at the point of transfer, since a small angle 
approximation, is used for the sine of y * 
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INSTRUCTION LISTING OF PROGRAM 2 


AV 

Di 

ai 

R$ 

ri 

R+ 

R+ 

4 

bi 

n 

s 

i 

b+ 

4 

aY 

5 
I 

b+ 

b \ 

s 


+ bT 

b+ c ■=■ 
B+ At 

a I — 

df A/ 

* ' 4 

bt /Z 

bi 

Bi /ff 

A-5- aZ 

b| 

AX aW 

At Bt 

- CT 

b$ ai 

BX ri 

c- R- 

b£ Rt 

BX Ri 

DX 


X CO 

A/ a- 

ai Di 

RO 

R+ A<> 

R- | 

-Ri bX 

Ri D$ 

rX CX 

Rt D| 

RS D- 

DX AO 

X c* 

4 c* 

/V D* 

V V 

aV 

A 0 

*0 

DO 


X 
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APPENDIX C 

PROGRAM 3 - LATERAL VELOCITY COMPONENT AND THE YAW AND 
PITCH ORIENTATIONS USING THE IN-PLANE VELOCITY 
COMPONENTS AND THE TOTAL VELOCITY COMPONENT AS INPUT 


ggggEDlNG 


PAGE BLAfci& 





APPENDIX C 

PROGRAM 3 - LATERAL VELOCITY COMPONENT ' AND THE YAW AND 
PITCH ORIENTATIONS USING T-HE IN-PLANE VELOCITY 
COMPONENTS AND THE TOTAL VELOCITY COMPONENT AS INPUT 

i „ ' 

The AV ' and the yay and pitch are calculated from AV. , AV„, and AV„. 

1 * • X A L 

This program requires a decimal wheel setting of 9 with V being the entry 
point of" the program. The order of the input and output quantities is 
as given . 



estrictions 

AV| If AV_j_ <_^AV^ + AV| , then the program will 
u. y^/ u v a negative number and return to the beginning. 


2. It must be noted that the AY 5 yaw and pitch output are absolute 

< JL 

values. This requires the user to place the proper sign on these values. 
Also, the yaw and pitch output takes on the following range of values: 
o 1* ±90, 0 < e < 90. As an example, if the yaw was output as 83°, but 
a retrograde component of velocity is desired, then the user knows the yaw 
desired is 97°.', A .similar situation holds true, for pitch. 






2k 


INSTRUCTION LISTING OF PROGRAM 3 


AV 

4 

bi 

A/ 

R- 

ai 

-bi 

ai 

B+- 

R* 

RO . 

AX . 

dS 

a/ 

Ri 

R+ 

C4 

c4 

A+ 

Ri 

R- 

c+ 

AW 

A+ 

R- 

R.4 

B$ 

a+ 

A+ 

R* 

•R4 

- t AX . 

d4 • , 

+ 

R4 

rX 

di 

Bi 

ai 

ri 

Ri 

B- 

AX 

ri 

D- . 

RS 

/z 

B- 

D4 

X 

DX 

AO 

a/ 

+ 

aO 

bi 

V 

ai 

A/ 

4 

S 

aZ 

d4 

ai 

/V 

4 

A/ 

Bi 

dX 

,C4 

bi 

bX 

A+ 

- 

di 

b| 

aO 

B+ 

A/ 

W 

s 

bi 

A/ > 

ai 

aV 

4 

AX 

B+ 

ri 

c* 

bi 

bi 

v 

d$ 

V 

AX 

c+ 

B+ ‘ 

X 


<4 

A/ 

A/ 

ai 


s 

b± 

B+ 

R4 
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APPENDIX D 

PROGRAM h - PITCH ORIENTATION BASH) ON THE TOTAL TELOCITY 
MAGNITUDE AND THE VERTICAL VELOCITY COMPONENT 



APPENDIX D 


PROGRAM 4 - PITCH ORIENTATION BASED ON THE TOTAL VELOCITY 

*2 

MAGNITUDE AND THE VERTICAL VELOCITY COMPONENT 


This program calculates pitch and the sine of the pitch as a function 

of AV. and AV_,. A decimal wheel setting of 9 is required for this program 
t a 

with V "being the entry point of the program. The order of the input and 
output is as given. 

Input ; AV_j_ , AV^ 

Output ; jsin fj, jp| 


Equations : 


sin 8 = 


AV, 


H= sin ‘ 1 (a^) 


Restrictions : 

1. AT Z < AV t 

2. The proper sign must "be placed on the pitch. 

3 . The range of values taken on "by pitch are : 0 <_ B £ 90 . As an 

example » if the pitch is output as 83°, hut a retrograde component of 
velocity is desired then the user knows that the pitch desired is 97°. 


HEEDING PAGE BLANK NOT FILMED 
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INSTRUCTION LISTING’ 'OF PROGRAM 4 


AV 

Bt 

Bt 

t 

Rt 

at 

<4 

A+ 

at 

R- 

RX 

S 

B+ 

rt 

R* 

RO 

t 

A/. 

Dt 

Rt 

R+ 

Bt 

. B+ 

+ 

Rt 

R- 

Ct 

A/ 

'a/ 

R- 

Rl 

aO' 

Bt 

at 

R* 

Rt 

at 

a/ 

dX 

Rt 

rX 

dt 

Bt 

- 

rt 

Rt 

Bt 

A/ 

A/ 

D- 

R5 

AX 

Bt 

at 

X 

DX 

B- 

A/ 

rt 

AO 

Bt 

a/ 

A+ 

Dt 

/0 

S 

at 

A+ 

X 

V 

4- 

dt 

A+ 

at 
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• APPENDIX E 

PROGRAM 5 - LATERAL VELOCITY COMPONENT AND ORIENTATION 
BASED ON IN-PLANE VELOCITY COMPONENT AND. THE 
TOTAL VELOCITY COMPONENT AS INPUT 



APPENDIX .E 

PROGRAM 5 - LATERAL VELOCITY' COMPONENT AND ORIENTATION 
BASED ON IN-PLANE ‘VELOCITY COMPONENT AND' THE 

t r t 

TOTAL VELOCITY COMPONENT AS INPUT 


The AVy, sin ifi, and are calculated from the sin 6, AV^., AV X * This 

program requires a decimal wheel setting of 9 with V .being the entry point 
of the program. The order of the input and output quantities is as given. 

Input : jsin b|, AV^. , AV^ 


Output : 



sin iM. 


1*1 


• Equations : 


AV Z = AVj. -sin 3 


m. 


:h V AV t - AT z - ?v x 


■Isin t[| = 


A V 


,'ipr. 


+ AV 2 


ij; = sin * f 


‘AV.. 


>K + av y 


Restrictions : 

1. AV 2 > AV 2 + AV 2 sin 2 £. If this condition is violated then the 

"C X b 

program will output a negative number and return to the beginning. 

. 2. The user should note that AV^, sin tp, and yaw are output as abso-i 

lute values , so this requires that he place the proper sign onto the 'com- 
puted value. The range of the value output for yaw -is 0 <_ ifj <_ 90. As an 
example, if the yaw was output as 83° , but a retrograde component of veloc 
ity is desired, then the user knows the yaw desired is 9T° 
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F t r !u, 'J2 e a l 0ma ' y ; ***** al,itu<)e ' 3nd perigee al,itude as a Unction of inertial velocity and 
inertial flight-path angle for various spacecraft altitudes. 
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Figure 2. -Continued. 
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(c) Altitude - 95 n. mi, 
Figure 2. - Continued. 
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Figure 2.- Continued. 
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Figure 2. -Continued. 
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Figure 2. - Continued. 
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(g) Altitude - 135 n. mi. 
Figure 2.- Continued. 
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(h) Altitude - 145 n. mi, 
Figure 2. -Continued. 
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Figure 2. -Continued. 
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Figure 2.- Continued. 
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Figure 2. -Continued. 
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Figure 2.- Continued. 
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Figure 2. -Continued. 



True anomaly, 0, deg 
Altitude of perigee, h^, n. mi. 

Altitude of apogee, h n. mi. 


, 100 - 


Perigee altitude., n. mi, 


j Apogee altitude, n. mi 


True anomaly , deg 




MBm 




aaBBfg 


-LO 


- 1.2 


-1.4 


- 1.6 


(n) Altitude 8 205 n. mi. 


24 840 24 880 24 920 24 960 25 000 25 040 25 080 25 120 25 160 25 200 25 240 25 280 25 320 25 360 25 400 25 440 

Inertial velocity, Vj, fos 


Figure 2.- Continued, 
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Figure 2.- Continued, 
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Figure 2. -Continued. 
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Figure 2. -Continued. 
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Figure 2. - Continued. 
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Figure 2.- Continued. 
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Figure 2. - Concluded. 
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( a ) v. = 24 800 - 25 900 fps; y. = .01 - 1.3 deg. 

Figure 3.- The Z component of velocity as a function of inertial velocity and inertial flight-path angle 
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(b) V = 24 800- 25 800 fps; y. = .1 - 1.3 deg. 
Figure 3.- Continued. 




Vj, fps 
25 300—3 


25 250H 


25 200 — 


25 150 — 


25 100H 


25 050-; 


25 000 -= 


24 950-; 


24 900 — 


24 850 — : 


24 800 



(c) V. = 24 800 - 25 300 fps; y. = .1 - 1.3 deg. 


Figure 3.- Continued. 
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Figure 3.- Concluded. 
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V., fps V x , fps Y jt deg 



(a) V. = 24 800 - 25 800 fps; y. = .1 - 1.3 deg. 

Figure 4.- The X component of velocity as a function of inertial velocity and inertial flight-path angle. 
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(b) V. = 24 800 - 25 300 fps; y. = .1 - 1.3 deg. 
Figure 4.- Continued. 
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(c) V. = 25 300 - 25 800 fps; y. = .1 - 1.3 deg 
Figure 4.- Continued. 
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(d) V. = 24 800 - 25 000 fps; y. = .01 - 1.3 deg. 
Figure 4.- Continued. 
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Ce) V. = 25 000 - 25 200 fps; y. = .01 - 1.3 deg. 
Figure 4.- Continued. 
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(f) V. = 25 200- 26 000 fps; y. = .01 - 1.3 deg 
Figure 4.- Concluded. 
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(a) AV total - 230 fps. 




Figure 5. - AVy as a function of pitch attitude and A V^; yaw attitude as a function of AV^ and AVy for various A V^, 
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(b) AV total - 330 fps. 
Figure 5. - Continued. 
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(c) AV total ■ 415 Ips. 
Figure 5.- Continued. 
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(d) A V total - 430 fps. 
Figure 5. - Continued, 



! AVy' , (PS 


1 



s 


(e) AV total ■ 475 fps. 
Figure 5. - Continued. 
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(f) AV total - 515 fps. 
Figure 5.- Continued. 
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(g) 6H total • 1100 fps. 
Figure 5. -Continued. 
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(h) AV total - 1200 fps. 
Figure 5. - Continued. 
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(i) AV total - 1300 fps. 
Figure 5. - Concluded. 


